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ABSTRACT. I present the motivation for studying nuclear collisions at ultrarelativistic energies which is to map 
the phase diagram of strongly interacting matter under very extreme conditions. The relevant experimental efforts 
are overviewed and problems in connecting observables with physics of the hot and dense matter are pointed out. 
I review the current highlights of the experimental programme: excitation functions of various quantities from 
the SPS which suggest a possibility of the onset of deconfinement, and the evidence of deconfined perfect fluid at 
RHIC. 



THE PHASE DIAGRAM OF QCD 

Phase diagram is graphical representation of the 
qualitative properties of bulk matter and of the condi- 
tions under which these properties change. Microscop- 
ically, the bulk properties result from the underlying 
interaction between constituents of the material. This 
also gives typical energy or temperature scales at which 
phase transitions would happen. 

All properties of materials we usually deal with 
are due to electromagnetic interaction. Let us take 
water as an example. Boiling and condensing is con- 
nected with rearranging bonds between the molecules, 
so intermolecular van der Waals interaction gives the 
typical energy scale, which is of order 10 -2 eV per 
molecule. By using typical molecular densities one can 
get an estimate for critical energy density and temper- 
ature. 

At somewhat higher energy the electrically neutral 
atoms break up into positively charged ions and elec- 
trons. This results in plasma phase. The typical ionisa- 
tion energy of an atom is of the order of electronvolt, 
thus about a factor 100 higher temperature than for 
boiling would be expected. 

We can ask if the strong interaction leads to ana- 
logic phase transitions. Strong interaction acts on en- 
tities with colour charge; the elementary participants 
are quarks and gluons. Similarly to electromagnetic in- 
teraction which acts also between electrically neutral 
atoms and molecules via van der Waals forces, there 
is strong interaction between colourless hadrons, e.g. 
nucleons within a nucleus. Nuclear matter is held to- 
gether by strong interaction. 

The typical energy scale of strong interaction be- 
tween nucleons in a nucleus shows up as the binding 
energy. In analogy to phase transition from liquid wa- 
ter to vapor known from everyday life, it turns out 
that there is a first-order "liquid-gas" phase transition 
in nuclear matter The typical temperature is of 

order of few MeV 1 . 

Is there a strong interaction analogy to electro- 
magnetic plasma phase with free charges? This ap- 
pears more problematic: while electric charges can exist 
alone in vacuum, colour charges cannot. They always 
must be bound together within a colourless object. 

iflere I use natural system of units where Boltzmann constant 
kg = 1 as well as Vi = c = 1. Temperature in SI units is obtained 
bv exDressine T first in Joules and then dividing: bv fcn. 



This property of quarks as carriers of colour charge is 
called confinement. Nevertheless, if the energy density 
is large enough, loosely speaking if the colour charges 
are packed together very densely, the quarks may not 
be anymore bound to hadrons but can freely roam over 
the whole region occupied by all quarks together 
Such a system is labeled "quark-gluon plasma" (QGP). 
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Fig. 1. Schematic view of the phase diagram of nuclear 
matter as a function of baryochemical potential and the 
temperature. Note that there are no scales put on the 
axes. 

In principle, the phase diagram is determined by 
the Lagrangian of underlying interaction. However, in 
case of the QCD Lagrangian it is technically too com- 
plicated to calculate from the first principles, yet. Its 
structure (Fig. is qualitatively and to certain level 
quantitatively known from effective theories and nu- 
merical simulations. The phase diagram is shown as 
characterised by the temperature T and baryochemical 
potential. Baryochemical potential is a measure of net 
baryon density. By changing temperature not only typ- 
ical momentum of quanta varies but also the densities 
of baryons and antibaryons. We must realise that the 
typical energy scales in this phase diagram are large, 
and relativistic effects of particle creation may appear. 
Thus the particle numbers are not conserved; what is 
conserved are quantum numbers like baryon number. 

At moderately high net baryon density there is the 
line of first-order liquid-gas phase transition. Farther to 
asymptotically high baryochemical potential while still 
at rather low temperature we find the regime of colour 
superconductivity 0]. This is an analogue to regular 
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weak attractive interaction between fermions on Fermi 
surface. Here, that role is played by certain channels of 
strong interaction. There are several different phases 
of colour superconductor, but we shall not go in such 
details here. Due to the low temperature required, such 
phase is not produced in nuclear collisions. It might be 
present in neutron stars, as they are rather cool: only 
some billion Kelvin. (If you think that this is hot, read 
below.) Still, there is no clear observational evidence 
for colour superconductor in a neutron star so far 

The regime I want to focus on is that of high tem- 
perature. We see in Fig. that at high temperature 
the plasma phase is expected. The equation of state 
at vanishing net baryon density can be calculated nu- 
merically in so-called lattice QCD approach (i . Such 
calculations indicate that at temperature 173±15 MeV 
the energy density changes as a function of tempera- 
ture dramatically but smoothly (Fig. [21 . In SI units 
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Fig. 2. The energy density normalised by T as a function 
of temperature from lattice QCD. Simulations with 2 and 
3 light quark flavours and with 2 light and one heavy 
flavour are presented. Labels show the energy densities 
reached at various accelerators. Figure from [Sj. 

this temperature is about 2 • 10 12 K, which is about 
100,000 times hotter than in the centre of the Sun! 
Two phase transitions happen here: confinement van- 
ishes and chiral symmetry is restored. The former is 
observed in the simulation as a change of quark-quark 
potential from linear rise to an asymptotically constant 
value at large distances; the latter as vanishing of vac- 
uum expectation value of chiral condensate (ijjip) ■ 

At large enough non-zero net baryon density first 
order phase transition appears. The critical line ends 
in a critical point of second-order phase transition. Ex- 
act position of this point is of fundamental interest, 
although so far there are only calculations based on 
effective theories or rather unreliable estimates from 
lattice QCD [7||E]. 

CREATING EXTREME CONDITIONS 

Have such extreme conditions ever been realised in 
nature? Can we create them artificially? The answer 
to both questions is positive. 

They existed in early Universe which evolved from 
a very hot and dense beginning through continuous 
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where Mp\ = ^Tic/Gm is the Planck mass and g* is the 
effective number of degrees of freedom (loosely speak- 
ing, how many different particle species contribute to 
the energy density). According to this relation, tem- 
perature around 150-200 MeV was reached at time 
10-100 microseconds after the Big Bang. This is the 
time when hadrons were born. 

The amazing thing is that we can re-create a small 
"early Universe" in the lab and study "Little Bangs" ! 
This is done with the help of ultrarelativistic nuclear 
collisions. By colliding large nuclei energy is converted 
into new quanta with large typical momentum. At the 
same time the system is large enough so that we can 
have bulk matter consisting of many interacting con- 
stituents. At sufficient collision energy the energy den- 
sity within the system will reach up to the regime of 
QGP. 

We have to realise that the colliding nuclei contain 
only baryons and no antibaryons, thus always matter 
with positive net baryon density is created. The higher 
the collision energy, the more secondary particles can 
be created. Among created particles there are as many 
baryons as antibaryons and so the ratio of baryons 
to antibaryons becomes smaller at higher collision en- 
ergies. Consequently, the baryochemical potential de- 
creases. 

Hence, nuclear collisions at various energies can 
map various regions of the phase diagram. Un- 
fortunately, such collisions create rapidly evolving 
systems — so-called fireballs — which makes it rather 
complicated to study the generated matter. Let us have 
a flash view of the evolution of fireball (Fig. 03). 
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Fig. 3. The space-time diagram of longitudinal coordinate 
and time of the evolution of fireball (a) without and (b) 
with the production of quark-gluon plasma. 



1. Colliding nuclei have relativistic energy, so before 
the collision they appear in the centre of mass 
system Lorentz contracted. The factor 7 spans 
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from about 4 (at lowest SPS energy) up to 100 
at RHIC and even more than 2000 at the LHC. 

2. If two nuclei collide at high energy, the incident 
nucleons and/or partons tend to continue their 
longitudinal motion while being just slowed down 
somewhat. After they have passed through each 
other, the energy they have lost converts into 
new quanta in the region between the receding 
original nucleons/partons. Just after the collision 
the energy density of this system is maximum 
and the matter could find itself in the deconfmed 
phase, if it is high enough. 

3. The longitudinal motion continues and thus the 
system expands longitudinally. Consequently, it 
cools down. Due to inner pressure also trans- 
verse expansion builds up which accelerates the 
cooling process. Hence, even if the matter was 
in QGP phase initially, eventually it changes to 
hadronic phase. This may happen via first-order 
phase transition in the part of the phase diagram 
at higher baryochemical potential or via rapid 
crossover. 

4. In any case, the fireball thus ends up in a 
hadronic phase where is still expands very fast. 

5. Eventually, the density of the system becomes so 
low that no more rescattering between its con- 
stituents appear and momenta of hadrons freeze- 
out. These are the momenta observed in detec- 
tors. It is important to realise that hadrons are 
observed and not quarks or gluons. 

This scenario runs very quickly. From the initial impact 
until the freeze-out it lasts about 10 fm/c or more 2 . 
Unfortunately, we are most interested in studying the 
hottest and densest phase of matter, but we observe 
hadrons escaping from the less excited hadronic gas. 
The main technical problem in this field is thus relating 
observable quantities with the properties of early hot 
state of matter. 

Before addressing this problem, let me list the exist- 
ing, recent, and forthcoming experimental facilities de- 
voted to this kind of experiments. For better overview, 
they are listed in Tabled Baseline in many observed 
effects has been provided by the Alternating Gradient 
Synchrotron (AGS) of Brookhaven National Labora- 
tory (BNL). This machine is currently used as preac- 
celerator for RHIC (to be mentioned later). Up until 
recently, experiments have been performed at CERN's 
Super Proton Synchrotron (SPS) which accelerated nu- 
clei up to the size of lead. These brought interesting 
results which could indicate the onset of quark-gluon 
plasma production in this energy region. Currently, 
world's strongest heavy ion accelerator is Relativistic 
Heavy Ion Collider (RHIC), where strongly interacting 
QGP behaving like a perfect fluid has been discovered. 
In less than a year from now the Large Hadron Col- 
lider (LHC) at CERN will be commissioned. Heavy 

2 The time unit fm/c is given by the time it takes for the light 

to pass the typical size of proton, which is 1 femtometre (called 

also fermi). This is the typical time scale in high energy physics. 
t„ ct ;4- ~~ — a~ ~u — . i- o in— 24 „ 



TAB. 1. Overview of accelerators used for 
ultrarelativistic heavy ion studies. I list the maximum 
energy per participating nucleon in the nucleon-nucleon 
centre of mass system and whether the setup is fix target 
experiments or it is a collider. 



machine 


lab 


/ „max 


target or collider 


recent 


AGS 
SPS 


BNL 
CERN 


4.8 GeV 
17 GeV 


target 
target 


present 


RHIC 


BNL 


200 GeV 


collider 


forthcomiE 


g 






LHC 
FAIR 


CERN 
GSI 


5500 GeV 
8 GeV 


collider 
target 



ion studies — in particular with the dedicated detector 
ALICE — are planned at this machine among others. In 
2014 completion of the Facility for Antiproton and Ion 
Research (FAIR) of the Gesellschaft fur Schwerionen- 
forschung (GSI) in Darmstadt is scheduled, which will 
deliver high luminosity ion beams with energies up to 
34 GeV per nucleon for fixed target experiments. 

In the following we shall mainly talk about results 
from experiments performed with lead or gold beams. 
There were experiments with smaller nuclei, as well 
(In, Cu, S, O, . . . ) in order to gain better systematics 
of various observables. 

FREEZE-OUT STATE: THE GROSS PIC- 
TURE 

The early hot phase of the collision is not observable 
directly. There are basically two approaches which can 
be used in order to learn about it. 

1. By studying momentum distributions of low p t 
hadrons one can extract information about the 
state of the fireball at the time of its breakup. 
This state can be reconstructed and one can de- 
duce the evolution of the fireball which could 
have lead to it. At least, this method can dis- 
criminate theoretical scenarios which lead to the 
wrong final state. 

2. Probes which are produced early and manage to 
escape the fireball carry direct information about 
the hot matter. These are e.g. particles which in- 
teract only electromagnetically, like photons and 
(di)leptons. Recently, hard jets have provided 
very clean probe of the matter, which we shall 
discuss below. 

Let us first have a look at what hadrons tell us about 
the final state. 

It is interesting to analyse their chemical composi- 
tion. This means that we study how many hadrons of 
different sorts are there among the observed particles. 
Interestingly, it turns out that the chemical compo- 
sition can be very reasonably described by chemical 



3 



equilibrium characterised by some value of tempera- 
ture and baryochemical potential. Moreover, the val- 
ues of temperature and chemical potential inferred for 
all collision energies lay on a well defined curve in the 
T-fx plane (Fig. EJ) fSJ. Results for RHIC and higher 
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Fig. 4. Temperatures and baryochemical potentials 
extracted from fits to chemical compositions of hadrons 
produced in nuclear collisions at various energies. Figure 
taken from ;9 . 

SPS energies seem to coincide with the region of rapid 
crossover from hadronic to partonic system and thus 
suggest that the abundances are fixed through the pro- 
cess of hadronisation. Chemical freeze-out at lower SPS 
energies, about 6 to 10 GeV per nucleon in the centre 
of mass system, might be close to the critical point 
of the phase diagram, though the precise position of 
this point is unknown, yet 0. This gives the first hint 
about regions of phase diagram which can be explored 
using different facilities: RHIC and LHC will look at 
hot matter at small baryochemical potential likely un- 
dergoing a smooth crossover to hadronic phase. On the 
other hand, SPS and the future FAIR can search for 
the line of first order phase transition and the critical 
point. 

From single-particle momentum spectra and mo- 
mentum correlations we know that the fireball ex- 
pands transversely rather fast ^JE]- In a transversely 
expanding fireball transverse momentum spectra are 
blueshifted due to Doppler effect. The inverse expo- 
nential slope in non-relativistic approximation is given 
as P2] 

= T + m(vt) 2 , where ; — oc exp [ — } . (2) 

p t dp t V T *J 

Thus from spectra of several identified species (with 
different masses) we can extract the average transverse 
collective expansion velocity (vt). In central collisions 
at RHIC it reaches up to 0.6c 

This figure is also supported by momentum correla- 
tions. They measure primarily the size of the emitting 



ergy in central Pb+Pb and Au+Au collisions the emit- 
ting fireball is about twice as big transversely as the 
original nucleus! The average momentum dependence 
of correlation functions is also consistent with strong 
transverse expansion Conclusion is that we 

do observe collectively behaving matter and not just a 
bunch of independent nucleon-nucleon collisions and 
that there is strong pressure in the early phase of the 
collision. 

PERFECT FLUID AT RHIC 

Interesting results have been observed in non- 
central collisions. Here, the two nuclei collide so that 
the reaction zone is not azimuthally symmetric but 
rather has elliptic shape with shorter size in direction 
of the impact parameter. This initial shape anisotropy 
is translated to an anisotropy in transverse momentum 
spectra: more particles are produced in the direction of 
impact parameter. This is quantitatively expressed by 
Fourier decomposition of the spectrum 

d 3 N _ 1 d 2 N 
p t dp t dy d<f> 2n p t dp t dy 

x (l + 2v 2 (pt,y)cos(<p-<j> r ) + ...) (3) 

(other terms of the decomposition vanish by symmetry 
at midrapidity, cj> r is the azimuthal angle of impact 
parameter). The second-order coefficient v 2 is called 
elliptic flow. 

Momentum anisotropy measured by the elliptic 
flow is very sensitive probe of conditions in the early 
hot phase. Data indicate that more transverse flow is 
produced in the direction of impact parameter than 
perpendicularly to it. This has natural explanation in 
hydrodynamics. Flow results from pressure gradients. 
The smaller size of the reaction zone in the direc- 
tion of impact parameter implies larger pressure gradi- 
ents and more flow. This mechanism of generating flow 
anisotropy vanishes with time as the fireball becomes 
azimuthally symmetric through expansion and so large 
flow anisotropy must be generated in the earliest phase. 

At RHIC, ideal hydrodynamics describes the trans- 
verse expansion and transverse momentum spectra (al- 
most) perfectly ^5]. (It only breaks for about 1% of all 
produced particles with highp t .) This has not been so 
at lower collision energies like SPS. Technically, hydro- 
dynamics consists of solving conservation laws and the 
equation of state. The simulation has been tuned on de- 
scription of central collisions, i.e., the prescription for 
determining initial conditions has been fixed there in 
such a way that spectra were reproduced. Then, with- 
out any other tuning elliptic flow in non-central colli- 
sions has been successfully predicted, as well ^B]- It 
should be stressed that the model was successful only 
if QGP phase was assumed during the early hot period 
of fireball evolution! 

Recall that data were reproduced with ideal hydro- 
dynamics, i.e. with vanishing viscosity. Different simu- 
lations with a kinetic model have shown that the mean 
free path must be extremely short, otherwise the sim- 
ulation would fall below the measured vo IT5I. Recall 
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Fig. 5. Azimuthal correlation function of high p t particles 
in pp, d+Au, and Au+Au collisions at i/s = 200 AGeV. 
Angle is oriented in direction of a high p t jet trigger 
particle. Figure taken from piij. 



The STAR collaboration at RHIC has studied az- 
imuthal correlations of high pt particles associated with 
jets (Fig. EJ) 21 . In proton-proton and deuteron-gold 
collisions they found indeed that jets were produced in 
pairs separated by 180° in azimuthal angle. However, 
in central Au+Au collisions the associated jet vanished 
(stars in Fig. EJ. Since momentum conservation still 
holds, the associated jet must have been produced but 
the leading particle has been completely slowed down 
when traversing the matter of the fireball. This is a sig- 
nature of large energy loss of leading jet parton which 
is only possible in a medium with high density of colour 
charge, such as QGP. 

Recently, more detailed studies showed that jets 
with very large transverse momentum poke out 
through the medium |22|. This opens possibility to 
study the energy loss of leading partons quantitatively 
and make conclusions about the density of early hot 
phase. 

THE ONSET OF DECONFINEMENT 



that in a kinetic approach shear viscosity grows with 
the mean free path because momentum can be trans- 
ferred to larger distances. It has also been shown that 
"viscous" corrections of the elliptic flow just due to 
modifications of locally thermal distribution are large 
and always decrease V2 E| • The conclusion is that the 
dimensionless 3 ratio of shear viscosity to entropy den- 
sity rj/s must be very small. Such a small viscosity 
is consistent with the picture of plasma with strongly 
interacting colour charges |18| . because weakly inter- 
acting gas would have long mean free path and larger 
viscosity. Calculations indicate that just above the crit- 
ical temperature the value of rj/s is as small 0.2 |19II2U| . 
For comparison, the corresponding value for water at 
normal conditions would be of order 4 10 3 . In fact, the 
strongly interacting quark-gluon plasma is the most 
perfect liquid known. 

DECONFINEMENT AT RHIC 

The previous section indicated that matter pro- 
duced at RHIC is deconfined. There is another spec- 
tacular demonstration of this: the jet quenching. 

Jets are well known in high-energy collisions of sim- 
pler systems. For example, in proton-proton collisions 
two partons from incoming protons can exchange large 
momentum and secondary particles are produced in 
two rather focused showers around the leading struck 
out partons. Due to momentum conservation, jets can 
never be produced alone. Most often, there is one as- 
sociated jet in the opposite direction, though in some 
cases three or more associated jets can appear and 
rarely even counterbalancing photon or lepton is emit- 
ted. 



3 This ratio is dimensionless in natural units. In SI, dimension 
of viscosity is [tj] = N ■ s ■ m -2 and that of entropy density is 
[s] = J ■ K" 1 ■ m" 3 , so [n/s] = K • s = [h/k B ]. 

4 I should say that at higher temperature rj/s of water can be 
as small as about 2, but even this is still much bigger than QGP. 



If QGP is produced at RHIC, what is the minimum 
collision energy at which deconfincment sets in? In or- 
der to answer this question one would study excitation 
functions, i.e. dependences of various quantities on the 
collision energy. 

One example is shown in Fig. H3 ,23]. The ratio of 
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Fig. 6. Excitation function of the ratio of multiplicities of 
positively charged kaons to pions. Curves present results 
of various model calculations. Data of the NA49 
Collaboration, figure taken from \23\. 

multiplicities of positively charged kaons to pions first 
rapidly grows as a function of collision energy. It peaks 
around y/s « 7 AGeV drops down a little and then it 
levels off. This interesting non-monotonic structure — 
which has been dubbed "the horn" — gives a first hint 
that something may be happening here. In addition 
to this, in the same energy region where the horn of 
Fig. El drops to the right one observes a plateau in the 
excitation function of the meanp t of kaons. Finally, the 
number of produced pions per participating nucleon 

cfiancps its slnnp as a function nf collision nnr-rov inst 
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at the same place as the horn |23| . 

These observations are successfully interpreted in 
framework of so-called Statistical Model of Early Stage 
24 which assumes the onset of deconfinement where 
the horn appears. However, if indeed QGP is produced 
here, then no hadronic model should be able to describe 
all data. Thus it is important to play Advocatus Diaboli 
and test hadronic models by comparing their predic- 
tions with data. Indeed, most of them fail in reproduc- 
ing these excitation functions EH1 123 ■ Recently, we 
constructed hadronic kinetic model |28| which repro- 
duced the horn of Fig.EI It remains to test this model 
on other pieces of data before it can be confirmed or 
ruled out. 

CONCLUSIONS 

Let me summarise my selection of current high- 
lights of ultrarelativistic heavy ion programme. 

1. QGP may be produced at lowest energies studied 
at the SPS. Careful investigations, also employ- 
ing the new FAIR machine, are necessary in order 
to confirm this. 

2. Jet quenching at RHIC implies that we have mat- 
ter which eats jets. The only known medium able 
to do this is quark-gluon plasma. 

3. Elliptic flow at RHIC is as large as it can pos- 
sibly be in perfect fluid and so the viscosity of 
QGP is tremendously small. This is consistent 
with the regime of strongly interacting QGP at 
temperatures just above T c . 

The Large Hadron Collider to be commissioned 
next year will shift the frontiers in points |21 and [21 
Due to higher incident energy, jets which probe the 
early hot matter will be produced copiously. At the 
larger initial temperature the medium could find itself 
in a not-so-strongly interacting regime as the strong 
coupling decreases for interactions at large momentum 
scale. This would lead to larger viscosity and it would 
be interesting to observe a corresponding decrease of 
elliptic flow. Such an observation would nicely extend 
our knowledge of the QCD phase diagram. 

REFERENCES 



7. M. Stephanov, Prog. Theor. Phys. Suppl. 153 
(2004) 139. 

8. Z. Fodor and S. D. Katz, Phys. Lett. B 534 (2002) 
87. 



20 



21 



27 



28 



9. J. Cleymans et ai, hep-ph/0607164 

10. B. Tomasik, U.A. Wiedemann, U. Heinz, Heavy 
Ion Physics 17 (2003) 105. 

11. J. Adams et ai. [STAR Collaboration], Phys. Rev. 
Lett. 92 (2004) 112301. 

12. T. Csorgo and B. Lorstad, Phys. Rev. C 54 (1996) 
1390. 

13. B. Tomasik and U.A. Wiedemann, Central and 
non- central HBT from AGS to RHIC, in Quark- 
Gluon Plasma 3, R.C. Hwa and X.-N. Wang eds., 
World Scientific, 2004, p. 715. 

14. M. A. Lisa, S. Pratt, R. Soltz and U. Wiedemann, 
Ann. Rev. Nucl. Part. Sci. 55 (2005) 357. 

15. U. W. Heinz and P. F. Kolb, Hydrodynamic de- 
scription of ultrarelativistic heavy ion collisions, 
in Quark-Gluon Plasma 3, R.C. Hwa and X.-N. 
Wang eds., World Scientific, 2004, p. 634. 

16. D. Molnar and M. Gyulassy, Nucl. Phys. A697 
(2002) 495; (Erratum ibid. A703 (2002) 893). 

17. D. Teaney, Phys. Rev. C 68 (2003) 034913. 

18. E. V. Shuryak and I. Zahed, Phys. Rev. C 70 

(2004) 021901. 

19. A. Nakamura and S. Sakai, Phys. Rev. Lett. 94 

(2005) 072305. 

A. Peshier and W. Cassing, Phys. Rev. Lett. 94 
(2005) 172301. 

J. Adams et ai [STAR Collaboration], Phys. Rev. 
Lett. 91 (2003) 072304. 

22. D. Magestro [STAR Collaboration], Nucl. Phys. A 
774 (2006) 573. 

23. C. Blume [NA49 Collaboration], Acta Phys. Hung. 
A 24 (2005) 31 |arXiv:nucl-ex /0405010|. 

24. M. Gazdzicki and M.I. Gorenstein, Acta Phys. Pol. 
B 30 (1999) 2705. 

25. H. Weber, E.L. Bratkovskaya, W. Cassing, H. 
Stocker, Phys. Rev. C 67 (2003) 014904. 

26. M. Wagner, A.B. Larionov and U. Mosel, Phys. 
Rev. C 71 (2005) 034910. 

Yu.B. Ivanov, V.N. Russikh, V.N. Toneev, Phys. 
Rev. C 73 (2006) 044904. 

B. Tomasik and E.E. Kolomeitsev, 
nucl-th/0512088| 



1. P.J. Siemens, Nature 305 (1983) 410. 

2. A.D. Panagiotou et ai, Phys. Rev. Lett. 52 (1984) 
496. 

3. J. C. Collins and M. J. Perry, Phys. Rev. Lett. 34 
(1975) 1353. 

4. M.G. Alford, K. Rajagopal, and F. Wilczek, Phys. 
Lett. B 422 (1997) 247. 

5. M.G. Alford, J. A. Bowers and K. Rajagopal, J. 
Phys. G 27 (2001) 541 [Lect. Notes Phys. 578 
(2001) 235]. 

6. F. Karsch and E. Laermann, Thermodynamics and 
in-medium hadron properties from lattice QCD, in 
Quark-Gluon Plasma 3, R.C. Hwa and X.-N. Wang 
eds., World Scientific, 2004, p. 1. 



